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Abstract 

Experimental data on the total cross section of e^e~ annihilation into hadrons are confronted 
with QCD and the operator product expansion using finite energy sum rules. Specifically, the 
power corrections in the operator product expansion, i.e. the vacuum condensates, of dimension 
d = 2, 4 and 6 are determined using recent isospin 7 = + 1 data sets. Reasonably stable results 
are obtained which are compatible within errors with values from r-decay. However, the rather 
large data uncertainties, together with the current value of the strong coupling constant, lead 
to very large errors in the condensates. It also appears that the separation into isovector and 
isoscalar pieces introduces additional uncertainties and errors. In contrast, the high precision 
r-decay data of the ALEPH collaboration in the vector channel allows for a more precise deter- 
mination of the condensates. This is in spite of QCD asymptotics not quite been reached at the 
end of the r spectrum. We point out that isospin violation is negligible in the integrated cross 
sections, unlike the case of individual channels. 
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1 Introduction 



In the framework of QCD sum rules [I]- [2], the quark and gluon propagator modifications 
due to confinement are parametrized in terms of vacuum expectation values of quark 
and gluon fields. These so-called vacuum condensates appear as power corrections to 
perturbative QCD in the operator product expansion (OPE) of current correlators at short 
distances. The QCD condensates are accessible from experiment through the spectral 
functions, i.e. the imaginary parts of the current correlators. Since from dimensional 
arguments there are no gauge invariant operators of dimension d = 2 built from quark 
and/or gluon fields in QCD, it has been customary to assume these power corrections to 
start at c? = 4. However, one cannot exclude a priori some dynamical mechanism capable 
of generating a term of (i = 2 in the OPE [3]. In a recent paper |1] it has been argued 
that there is no evidence for such a term in lattice QCD, contrary to earlier claims. There 
is also no definite evidence from r-decay as shown in recent sum rule analyses |S]-[H] 
based on the final ALEPH data [9]. These determinations have also found values for the 
dimension d = 4 and d = Q condensates consistent with expectations. 

Given the impact of these condensates on QCD sum rule applications it is important to 
attempt a determination based on experimental information independent of r-decay. His- 
torically, determinations based on e~^e~ annihilation data have preceded r-decay analyses 
[To] . Given the current availability of more accurate data from e~^e~ annihilation, and 
the sizable increase in the value of the strong coupling over the years, it is imperative to 
perform a fully updated determination of the condensates to confront with results from 
r-decay. An additional motivation for such a project is related to the current discrepancy 
in the hadronic vacuum polarization value of the g — 2 factor of the muon obtained from 
r-decay and from e'^e~ annihilation data jll], to wit. The Standard Model (SM) pre- 
dictions of a^, the anomalous magnetic moment of the muon, are limited in precision by 
contributions from the hadronic vacuum polarization. The dominant terms are usually 
calculated via a weighted dispersion relation involving experimental cross section data 
on e'^e~ annihilation into hadrons at low energies, and perturbative QCD (PQCD) at 
high energies. The required integration kernel emphasizes low photon virtualities. Dur- 
ing the last decade the SM prediction of has improved continuously as more precise 
experimental data became available [12]. Nonetheless, a discrepancy between the mea- 
sured and the predicted values remains and it is tempting to trace this discrepancy back 
to new physics beyond the SM. Unfortunately, the experimental hadronic data entering 
the analysis originate from different experimental groups and often errors do not overlap. 
In addition, most of the modern experiments are dominated by systematic uncertainties 
which are difficult to estimate. Consequently, one should test the data base to the great- 
est possible extent against the underlying theory, i.e. QCD, before accepting the drastic 
conclusion that the SM has failed to describe the data. 

In the present paper we attempt to investigate the status of experimental e'^e~ anni- 
hilation data by confronting them with the OPE using QCD sum rules. We study the 
first three moments Mn^sq) {N = 0, 1, 2) of R{s), the ratio of the cross sections for 
hadron to muon production in e'^e~ annihilation, at a finite upper limit Sq of the integra- 
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tion over the square of the center-of-mass energy, sq. These moments are related to the 
QCD condensates of dimension d = 2N + 2. The extracted values of these condensates 
are very sensitive to experimental and theoretical errors as they result from evaluating a 
small difference between two large numbers, i.e. the difference between integrated data 
and PQCD. The sum rules therefore constitute a sensitive test of the consistency of the 
data with general expectations from QCD. In parallel with the determination of a poten- 
tial d = 2 term in the OPE, we also consider power corrections of dimension d = 4 and 
d = 6. This can help establish whether there is evidence for a dynamical d = 2 power 
correction, and whether the results for the higher dimensional condensates are consistent 
with expectations. There is also independent information on the d = A gluon condensate 
from charmonium analyses, as well as theoretical estimates of the d = 6 four-quark con- 
densate which can be used to check for consistency. We limit the analysis to dimension 
d = 6, as higher dimensional condensates are affected by such large uncertainties that no 
meaningful value can be obtained at present. 

In the evaluation of the moments, the upper limit of the sum rule integral must be taken 
large enough for PQCD to be applicable. It should be recalled in this connection that the 
QCD asymptotic regime is reached in the time-like region only at energies much higher 
than in the space-like region. For instance, this is seen in r-decay data of the ALEPH 
collaboration [5], where the V — A hadronic spectral function does not agree with the 
vanishing result of PQCD even at sq ~ fn^. In contrast, deep inelastic electro n-nucleon 
scattering is well known to exhibit precocious scaling. 

In connection with the determination of the muon magnetic moment anomaly, a^, only 
the sum of the isospin 1 = 1 and / = spectral functions is needed. However, QCD 
also makes definite predictions for the moments of the separate isovector and isoscalar 
components of the spectral functions. The corresponding moments are related by a trivial 
overall factor determined by the quark charges. Therefore, from the theoretical point 
of view it is interesting to study both the isovector and the isoscalar channels of R{s) 
separately. Of specific interest is the 1 = 1 channel since independent data are available 
in this case from hadronic r-decay. For example, it was suggested in [13] to use the r- 
spectral functions in the calculation of the muon anomaly. To be precise, the authors of 
[T3] replaced the r — ?■ z^r + tt^tt^, 27i^7i^7i'^ and tt^Stt^ spectral functions measured in r- 
decay with the poorly known neutral 27r and 4tt spectral functions from e~^e~ annihilation. 
For these exclusive final states, sizable model-dependent isospin-breaking effects need to 
be incorporated. This is in contrast to the case of moments, as opposed to individual 
channels, where isospin breaking is negligible for large enough sq. But if the aim is to 
determine QCD condensates for definite isospin, then the necessary separation of e'^e~ 
data into final states with / = or / = 1 would require a detailed understanding of 
isospin violation. This separation, though, would not even be possible in the presence 
of mixing. This mixing has been known for quite some time to be sizable in the p-u 
system. In addition, it was shown recently that p-7 mixing could, perhaps, explain the 
well-known discrepancy between e~^e~- and r-decay-based evaluations of [H]. However, 
for a determination of QCD condensates one does not need such a detailed, and generally 
model dependent analysis. In fact, isospin is almost exactly conserved in integrated cross 
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sections or moments, provided the upper limit of the integration over squared energies 
energies, Sq, is taken large enough. This is because isospin-breaking terms are proportional 
to m'^/sQ where rriq refers to the mass of the light quarks. Therefore one can restrict the 
analysis of moments and condensates to the total (/ = 1 plus / = 0) e+e~ cross section. 

A similar analysis based on the known 1 = 1 spectral functions obtained from r-decay 
data has already been performed [5]-[8]. However, there the upper limit of integration sq 
is restricted by the r mass to values of about sq — 3 GeV^ and it is not clear whether this 
value is large enough for PQCD to be applicable. In spite of this, the condensates derived 
from r-decay have been shown previously to be consistent with QCD expectations and 
we will therefore focus in the following on the case of e~^e~ annihilation. 

The paper is organized as follows. In Section 2 we present a short summary of the 
theoretical background and in Section 3 we describe shortly the data used in our analysis, 
with a more detailed discussion deferred to the Appendix. Section 4 contains our results, 
their interpretation, and our conclusions. 



2 QCD Sum Rules 

We begin by considering the electromagnetic current correlator 

n™(g2) = .|rf^xe^''^(0|T(J™(^)jr(0))|0) 

= (-gV + g.g.) n™(g^), (i) 

where for three flavours 

Jj^'^ix) = lu{x)Yu{x) - ld{x)Yd{x) - ls{^)Ys{x) . (2) 

QCD is flavour blind, and if isospin invariance is exact, it is convenient to define a QCD 
current correlator using any of the quark currents qi'J^qi with flavour i. This leads to 

Imn™(g2) = J]g^Imnvv(g^), (3) 

1=1 

where Qi is the charge of the quark i = u, d, s, . . ., = 2/3 foi nj = 3 flavours, and 

Imllvvl?^) is the QCD correlator of vector currents of flavour i. The electromagnetic 
spectral function, Imn^'^(s), with s the square energy, is accessible experimentally from 
data on e~^e~ annihilation to hadrons as follows. The ratio R{s) is defined as 

s crTOT(e+e" hadrons) 
= . , _ , , „^ , (4) 



with 



a{e^e /i^/i = — , (5) 

3s 
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and OEM = e^/in. In QCD the ratio R and the electromagnetic spectral function are 
related through 



Uf 

R{s) = 12 vr Im n™(s) = 3 5^ (l + ^ + . . .) . (6) 



A singlet contribution proportional to (^,. Qi)'^ arises at order 0{a^) and vanishes if one 
sums over three flavours. In particular, for the two-pion final state, dominated by the 
p-resonance, the relation between R and the pion form factor, Fi^\s), is 

Re^e^^.^.-is) = 7 - \Fi'\s)\' . (7) 



On the QCD side, the vector correlator is assumed to satisfy the OPE extended beyond 
perturbation theory. Non-perturbative modifications due to confinement are parametrized 
in terms of vacuum condensates, i.e., 

oo ^ 

87rnvv(g^) = E 77^^2^(^'""') (0|(^2^(/^')|0) , (8) 

where = —q^ > is large, /i is a renormalization scale, and the first term with 
N = stands for the PQCD contribution. The Wilson coefficients C2N, calculable in 
perturbation theory, contain the short distance information while the vacuum condensates 
effectively parametrize the long distance dynamics. These condensates are organized 
according to their dimension, with the leading ones being the product of the quark mass 
and condensate of dimension d = 4, mq{0\qq\0) , and the gluon condensate also oi d = 4, 
{0\^G^v G^IO). This gives the result 

2 

^4(04) = ^M-^n+4 7r2 ^m^(^uu)+m,{dd)+ms{ss)) , (9) 

where is the running strong coupling, and in the sequel (0|C2Ar|0) = {O2N) is to be 
understood. This condensate is renormalization group invariant to all orders in PQCD. 
Next, at dimension d = 6 there enters the four-quark condensate 

Ce{Oe) = -AiT^as{u-f^^5t''u-d-f^^5t''dY) 

-^7r3(a, (n7^rn + (J7^t"d) qif^t^'q) , (10) 

q=u,a,s 

which has a mild dependence on the renormalization scale. A once popular approximation 
is that of vacuum saturation [2] which gives 

8Qfi 

Cemivs = -^^'c^s\{qq)\' ^ -0.025 GeV^ . (11) 

This value serves only as an order of magnitude reference, as there is no reliable way of 
estimating corrections to this approximation. In our analysis we do not need to invoke 
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vacuum saturation as we will be determining the complete C^{Oq), independently of any 
theoretical assumption. While a dimension d = 2 contribution is present in the PQCD 
term of the OPE for the vector correlator, it is numerically negligible as it is proportional 
to + m^. 

The PQCD vector spectral function is well-known to five- loop order [I5]-[in], i.e. 



Svrlm nvvl-s) 



1 + a, + f F3 + 

vr' „o _ 5 



/33 



(12) 



where a, = a,(/i2)/7r, = ln(gV/i2), = 49.076 [20], F3 = 1.9857 - 0.1153n/, F4 



18.2427 



3 ^ 2 



- 4.2158n/ + 0.0862^2, = -1(11 - |n/), (32 



1(102 - fnj) 



= — 32(^^ ~ 'tI^"'/ + W'^/)' ^^"^ running coupling to four- loop order is [2T 



TT 



TT 



TT 



/3i 



+ 



TT 



(In^L-lnL-l) 

Pi 



/3i 



«s (so) 



TT 



^ f ln=^L - -In^L - 21nL ^ 

Pi 



(13) 



with 



ai^^(so) 



TT /3iL 

where L = ln(so/A|^) defines the standard MS scale Aj^, and the constant 63 is 



(14) 



1 

44 



149753 /1078361 6508, , 

+ 3564(3 - + ri^ 



6 

/50065 6472 
V 162 ^ "sT 



C3 + 



162 

1093 
729" 



27 



-n. 



(15) 



with C3 = 1.202. 
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Invoking Cauchy's theorem in the complex squared energy s-plane allows to relate the 
experimentally measured hadronic spectral function with that from QCD (quark-hadron 
duality), leading to the finite energy sum rules (FESR) 

i-f C2N+2 {O2N+2) = Stt' r ds - Imn°^T^(s) - s^+'M2N+2{so) , (16) 

Jo 

where the dimensionless PQCD moments M2n+2{so) are given by 

N 

-Imn(s)PQ^°. (17) 

It has been assumed that sq is large enough so that the replacement n(s) — 11^^^° (s) 
is justified on the circle of radius |so| in the complex s-plane. It must be mentioned 
that there is no mixing of operators of different dimension up to second-loop order in the 
coupling [22], hence each FESR involves only one condensate. We shall neglect radiative 
corrections to the condensates and proceed to determine the condensates of dimension 
d = 2 to d = 6 using the FESR, Eq. (fT6|l . This approximation is justified a-posteriori, 
given the resulting large errors in the condensates. 

The contour integral in the complex s-plane will be evaluated in the framework of Fixed 
Order Perturbation Theory (FOPT), as well as in Contour Improved Perturbation Theory 
(CIPT). In FOPT the coupling and the quark So are considered fixed 

(constant) so that only logarithmic terms contribute to the integral. The renormalization 
group (RG) summation of leading logs is only carried out after the contour integration 
by setting /i^ = — sq. In this case the integration can be done analytically. In the case of 
CIPT the strong coupling and the quark masses are running and the RG is implemented 
before integrating. The RG equation for the running coupling and quark masses is solved 
numerically at each point on the circle of radius sq. The value of the strong coupling to be 
used here is asiM^) = 0.321 ± 0.015, corresponding to A^g = 341 ± 24 MeV in FOPT, or 
as{M^) = 0.344 ± 0.014 corresponding to A^g = 382 ± 24 MeV in CIPT |23]. If one were 
to compare results for the condensates obtained from e+e~ data as described above with 
corresponding results from r-decay, one would need to change the overall normalization 
of the QCD correlator. 



/ N ^ 9 1 / ds s 

M27V+2 So = -Stt'— (b 

27rz Jc(u\) So _So_ 



n ds 


s 


Jo So 


So 



3 Hadronic Data 

The most rigorous approach to making a collection of exclusive data is to combine all 
available data for a given final state [2l]-[5T]. Their respective statistical and systematic 
uncertainties should be used to define the weight in the averaging procedure. In fact, 
this is important for precision determinations, such as, e.g., the g — 2 factor of the muon. 
However, one can be less rigorous for the purpose of the present analysis. The philosophy 
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Number 


e e — > 




'7" w 1 n2/T^^"\7'4\ 

Xi X lu (^v^ev ) 


1 




1.51 (4)(9) 


0.99 (3)(6) 


2 


777 


0.30 (4)(2) 


0.29 (5)(2) 


3 




133.1 (25)(8) 


95.6(27)(7) 


4 




30.1 (3)(14) 


41.6(8)(16) 


5 


2(7r+7r-) 


60.5 (3)(30) 


156.6 (10)(78) 





7T 7T Z[7T ) 




ZUD.D (^loj^lol ) 


7 


27r+27r-7rO 


7.1(2)(7) 


20.8 (6)(20) 


8 


3(7r+7r-) 


2.59 (7)(2) 


8.83 (30)(60) 


9 


2(7r+7r-7r") 


11.6(3)(12) 


39.2(10)(41) 


10 


77(7r+7r^) 


6.9(3)(6) 


18.9(8)(16) 


11 


rju! 


3.1(1)(3) 


9.0(3)(9) 


12 


r]{2TT+2TT~) 


1.0(1)(1) 


3.6(4)(3) 


13 




2.22 (4)(6) 


4.5(1)(1) 


14 




25.2 (3)(6) 


34.0 (7)(9) 


15 




12.2 (2)(5) 


14.2 (5)(7) 


16 




3.6(2)(3) 


11.5 (7)(8) 


1 7 






3 ^^)^'2) 


18 


nn 


0.9(3)(1) 


3.5(13)(5) 


19 


inclusive 


55.0 (4j(2j 


231.6 (30j(9) 


20* 


TT+TT- (CHPT) 


0.01(0)(0) 


0.001(0)(0) 


21* 


7r^7r^37r'^ 


3.5(1)(12) 


10.4(3)(34) 


22* 


TT+TT- (47rO) 


1.9(1)(6) 


6.2(2)(21) 


23* 


WTT+TT", w2tt'^(ci; — !> TT*^7) 


0.56(2)(18) 


1.6(1)(5) 


24* 


a;(non-3TT, TT7, 777) 


0.11 (1)(4) 


0.07 (1)(2) 


25* 


(j>{TiOTi-KK, 37r, TT7, 777) 


0.04 (0)(1) 


0.04 (0)(1) 


26* 


77TT+Tr-(2TT°) 


1.0(1)(3) 


3.6 (4) (12) 


27* 


KKtt 


16.6(4)(57) 


47(12)(16) 


28* 


KK2tt 


36.1 (9)(64) 


124 (3) (22) 


29* 


KK3tt 


1.7(1)(3) 


6.1(5)(10) 




Total 


499 (3)(13) 


1105 (5)(38) 



Table 1: The contributions to X^i = ^ Jq° s"'R{s) ds, using sq = 4.5 GeV^. The first error is 
statistical, while the second is systematic. For the total uncertainty, we added the uncertainties 
of all the final states in quadrature. Stars (*) indicate that the cross section had to be estimated, 
following the prescriptions explained in Appendix 
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^0 (GeV^) 



Figure 1: C2 {O2) calculated in FOPT using as{M'^) = 0.321 ± 0.015, corresponding to A^"^^ = 341 
MeV. The smaller uncertainties are obtained assuming no correlations between experiments, while the 
larger ones assume 100 % correlations for data obtained using the same experimental facility. 

is to make use of only the most recent data for a given exclusive hadronic final state. 
If the most recent data do not cover the entire energy range for that final state, then 
older data are used to cover this region. In this way there are no overlapping data 
sets, and then one does not need to be concerned about how precisely different, often 
inconsistent, data are combined. We will make use of exclusive data up to a/s = 2 GeV, 
and above this energy we use the BES inclusive data These data are consistent 

with PQCD, hence one observes a plateau for the condensates above ^/sq = 2 GeV. Since, 
in addition, their systematic errors are small (3 %) we observe that the resulting errors 
for Cn {On) do not change any more above this energy. The treatment of the statistical 
uncertainties is standard. We assume no correlations either within individual experiments 
or between different experiments. The exception to this is the treatment of the tt+tt^ 
data from BaBar [23]. In this case we are given a full correlation matrix. To simplify 
the calculation, we simply assumed 100 % correlation for the statistical uncertainty, which 
gives an upper bound on the uncertainty^] For the systematic errors, we assume that they 
are completely correlated within a data set. It is reasonable to assume that there will also 
be correlations between different final-state data sets obtained by the same experimental 
facility, e.g. BaBar, CMD2 or SND, since each experiment would use the same procedure 
when applying radiative corrections or for the determination of the luminosity. As we 
do not have any information about these correlations, we will quote our final uncertainty 
both for the worst-case scenario assuming a 100 % correlation of data determined at the 
same facility, as well as for the best-case scenario of no correlations. The true uncertainty 
can be taken to lie between these two cases. For more details see the Appendix where we 
also present the list of all final states used in our analysis in Table HJ We have verified our 

^We also checked that this assumption does not have a significant impact on the uncertainty in the 
final results for the condensates. 
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^0 (GeV^) 



Figure 2: C4 {Oa) calculated in FOPT using as{M^) = 0.321 ± 0.015, corresponding to A^^"^^ = 341 
MeV. The smaller uncertainties are obtained assuming no correlations between experiments, while the 
larger ones assume 100 % correlations for data obtained using the same experimental facility. 
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Figure 3: Cg (Oe) calculated in FOPT using a,(M2) = 0.321 ± 0.015, corresponding to A^"^^ = 341 
MeV. The smaller uncertainties are obtained assuming no correlations between experiments, while the 
larger ones assume 100 % correlations for data obtained using the same experimental facility. 

data base against the one used in [TT] for the determination of the muon anomaly, a^. 
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4 Results and Conclusions 



Results for the condensates of the total / = plus 1 = 1 current correlator determined 
from e^e~ data are presented in Figs. [TH From standard QCD and the OPE one would 
expect the dimension d = 2 condensate to be zero. Fig. [1] shows that this is the case 
within (large) errors, although the central value is negative. It is also seen from this 
figure that the sum rule starts to be saturated only at about sq ~ 3.4 GeV^. For the 
condensate of dimension (1 = 4 shown in Fig. [2] we observe saturation of the sum rule at a 
slightly higher sq than in the case of dimension d = 2. The central value of the dimension 
d = 4 condensate is seen to be mostly positive, as it should be since it is related to the 
vacuum energy. Figure [3] shows the result for the dimension d = 6 condensate. Since the 
errors are so large, no meaningful conclusion can be drawn in this case, other than that 
the central value is negative, in agreement with the vacuum saturation approximation, 
Eq. fllip . The results for the condensates are not very sensitive to the value of used in 
the sum rules. Their values at sq = 4.5 GeV^, assuming uncorrelated data errors are 



C2 {O2) 



C4 (O4) 



Ce (a) 



-0.13±0.13|exp±0.03| 
-0.14±0.13|exp±0.03| 

+0.16±0.38|exp±0.04| 
+0.10±0.38|exp±0.04| 

-0.40±1.30|exp±0.10| 
-0.40±1.30|exp±0.10| 



GeV^I 



FOPT 



GeV^lciPT 



GeV^I 
GeV^I 



GeV^ 
GeV*^ 



FOPT 
CIPT 

FOPT 
CIPT 



(19) 



(20) 



Assuming instead maximally correlated data errors the experimental uncertainties above 
increase to ±0.18|exp, ±0.52|exp, and ±1.70|exp, for C2 (C2), C4 {O4), and Cq (Oq), re- 
spectively. These results are consistent with expectations and with the values extracted 
from r-decay [S]-[H], albeit within very large errors. The uncertainties in all the results are 
clearly dominated by the experimental errors. One consequence of this is that e.g. in the 
d = 2 case it is not possible to distinguish among various theoretical sources of such a term 
in the OPE. This is also the case with the result from the r decay data [5]- [6], although it 
has a smaller uncertainty. In any case, the impact of the PQCD input for a given integra- 
tion method (FOPT or CIPT) is such that the higher the value of ctg the lower the curve 
in FigJTJ A similar trend is found from additional higher order PQCD terms, and this 
behaviour is common to all condensates. For the gluon condensate the value determined 
from charmonium data is C^iO^) = 0.06 ± 0.04 GeV^ [2], [52]-[5l], where the error is our 
conservative estimate. The result from r-decay is [6] €4(04) = 0.07 ± 0.02GeV'^ for a 
similar input value of a^. Concerning the d = 6 result, given the large uncertainty in its 
value no meaningful comparison can be made with e.g. the vacuum saturation predic- 
tion, Eq. lITT]) . It should be emphasized that the consistency between data and QCD is 
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Number 


e^e — > 


Reference 


Data Range 
Vs (CreVj 


Vacuum Pol. 


1 




CMD-2 (2005) [H] 


0.6-1.31 


Dressed 


2 




CMD-2 (2005) [21] 


0.6-1.38 


Dressed 


3 




BaBar (2009) [25] 


0.31 - 2.95 


Bare 


4 




BaBar (2004) [55] 


1.06-2.99 


Dressed 






CMD-2 (2006) [27] 


1.01 - 1.06 


Dressed 






SND (2002) [58] 


0.98-1.01 


Dressed 






SND (2003) [ig 


0.66-0.98 


Dressed 


5 


2(7r+7r-) 


BaBar (2005) [30] 


0.62 - 4.45 


Bare 


6 


7r+7r-2(7rO) 


BaBar (2010) [3T] 


0.76-3.31 


Dressed 






SND (2009) [32] 


0.66-0.94 


Dressed 


7 


27r+27r-7rO 


BaBar (2007) [33] 


1.03-2.98 


Dressed 


8 


3(7r+7r-) 


BaBar (2006) [34] 


1.3-4.5 


Dressed 


9 


2(7r+7r-7rO) 


BaBar (2006) [M] 


1.3-4.5 


Dressed 


10 




BaBar (2007) [53] 


1.03-2.98 


Dressed 


11 


rju) 


BaBar (2006) [33] 


1.25 - 3.25 


Dressed 


12 


r]{2TT+2TT-) 


BaBar (2007) [33] 


1.31 - 2.89 


Dressed 


13 


W7r°(a; -> 7r°7) 


CMD-2 (2003) [35] 


0.92-1.38 


Bare 


14 




CMD-2 (2008) [35] 


1.011 - 1.034 


Dressed 






Qivfn /'o^^7^ 

oiNU I^ZUU/j [oY j 


10/1 1 
1.U4 — 1 .00 


Dressed 






DM2 (1988) [3H] 


1.38-2.40 


? 


15 




SND (2001) [39] 


1.01 - 1.06 


Dressed 






SND (2006) [40] 


1.04- 1.38 


Dressed 






DM1 (1981) [H] 


1.4-2.1 


? 


16 


r](j) 


BaBar (2008) [42] 


1.57-3.45 


Dressed 


17 


PP 


BaBar (2006) [13] 


1.88-4.2 


Dressed 


18 


nn 


Fenice (1998) [H] 


1.9-2.44 


Dressed 


19 


Inclusive 


BBS (2002) [45] 


2-5 


Bare 






BBS (2009) [ig 


2.6-3.65 


Bare 


on 


TT^TT (XPl) 


* 


U.i4 — U.oi 




21 


7r'''7r~37r° 


* 


1.03-2.98 




22 


7r+7r-(47rO) 


* 


1.3-4.5 




23 


tt'7r"*'7r~. w27r'^(aj — > 7r'^7) 


* 


1.15-2.53 




24 


tt'(non-37r, 7r7, 777) 


* 


0.7-0.8 




25 


(j}{TLOTi-KK, 37r, 7r7, 777) 




1.01 - 1.03 




26 


T77r+7r-(27r°) 


* 


1.3-2.9 




27 




* 


1.3-4.7 




28 


if;^27r 


* 


1.4-4.3 




29 


ifX37r 


* 


1.6-4.5 





Table 2: Data used in evaluating the condensates. Stars (*) indicate that the final state has 
not been measured and was estimated as explained in the text. The treatment of the vacuum 
polarization correction is described in the text. Question marks indicate that it is not known 
which corrections have been applied by the experimental collaboration, and no corrections were 
applied by us in these cases. 
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in no way obvious. The high sensitivity of the FESR can be demonstrated if one were 
to determine condensates from the individual / = and 1=1 channels separately. For 
example, for the non-strange 1 = 1 final states listed in lines 2, 3, 5, 6, 8-10, 12, 22, and 
26 of Tables [T] and [21 the dimension d = 2 and d = 4 condensates would turn out to 
be unacceptably large, i.e. more than 3a away from theoretical expectations. Moreover, 
with this subset of data there would no longer be a clear plateau in the so-dependence , 
as the one in Figs. [H [2] and [31 We believe that these unacceptable results indicate that 
the separation of the data into 1 = 1 and 1 = components requires a fairly detailed 
understanding of isospin symmetry breaking. If isospin symmetry is broken, states with 
different isospin quantum numbers can mix and unknown interference contributions can 
spoil a naive decomposition into a sum of isovector plus isoscalar parts. 

To summarize, with the help of FESR we have confronted QCD and the OPE with 
experimental data for R{s). In particular we have considered the zeroth, first and second 
moment of R{s). This comparison of theory with experiment can simultaneously help to 
give support to the theoretical framework, i.e. the basic principles of QCD, the OPE and 
analyticity, and as a test of the quality and completeness of the data. We found that recent 
e~^e~ annihilation data are consistent with QCD and the OPE within very large errors, 
even though there exists some indication that not all data are yet correctly accounted 
for. Any future improvement on the accuracy of vacuum condensate determinations from 
e~^e~ annihilation will require a considerable reduction in the experimental uncertainties. 

A Hadronic Data 

In this Appendix we give details of the data sets used in our analysis. Table [2] lists all 
relevant final states, their corresponding energy ranges and the corresponding references. 
Some of the exclusive channels are as yet unmeasured. We give here the estimated values of 
these missing channels. These estimates are mostly based on isospin arguments, primarily 
as described in Refs. [TT] and |17j, where further details can be found. It should be 
mentioned that the Pais isospin-class analysis [TT] can be employed only after removing 
the ?7- contribution from some of the final states. This also avoids double counting in 
the evaluation of the FESR integrals. We have only quoted the results. A 33 % model 
error has been assigned to all of these data, unless stated otherwise. The numbers in the 
following list refer to the final states listed in Table [H 

20. There are no data for the 7t~^tt~ final state at low energies, so the prediction of 
chiral perturbation theory (CHPT) has been used. The cross section is expressed in 
terms of the vr form factor, cr(7r"'"7r~) = vra^/^Q |F°p/(3s), where /3o = (1 — 4m^/s)"^/^. 
The form factor is taken of the form = 1 + ^{r'^).,^^ + cis'^ + C2S^ + 0{s^). The 
constants are found by fitting the form factor to space-like data, with the result 
(r^)^ = (0.439 ± 0.008) fm^ a = (6.8 ± 1.9) GeV"^ and C2 = (-0.7 ± 6.8) GeV"® 

m- 
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21. A prediction for the missing channel cr(7r+7r 37r°)^_.oxci. = |cr(27r+27r n^)r)~cxc\. can 
be obtained after removing the ?7-contribution: cr(27r+27r~7r°)^_cxci. = o'(27r+27r^7r°) — 
a{7]7r+n~) x B{r] TT+n~n^), where B{r] 7r+7r-7r°) = 0.2274 ± 0.0028 and where 
a(27r+27r~7r°) is obtained from BaBar (2007) [SB]. 

22. cr(7r+7r-47r°) = 0.0625 a(37r+37r-) + 0.145 a(27r+27r-27r°)^_cxci from isospin symme- 
try, where the ?7-subtracted cross section (T(27r"'"27r~27r°)^_exci is determined using 
cT(27r+27r-27r°)^_exci. = a{2Tr+2Tr-2n°)- a {r]uj)xB{r] n+7r-n°)x B{uj tt+tt-tt^), 
where B{lu tt+tt~tt^) = 0.892±0.007, a(27r+27r-27rO) was measured by BaBar [M], 
and (T{riu) was taken from [M] . 

23. One has cT(a;7r"^7r~) from BaBar (2007) [33]. However, the dominant three-pion 
decay of the u already appears in the five-pion final state. Thus we calculate the 
contribution for uj — )■ -7r°7 by using a^uir'^'iT^ 7r'^7r"7r°7) = a{uj'iT'^7T~) x B{uj — > 
71^^) where B{uj — t- tt^^) = 0.0828 ± 0.0028. Then from isospin arguments it follows 
that cr((jj27r° — 37r°7) = 0.5 a{ujiT^iT^ — t- 7r+7r~7r'^7). 

24. There are cross section data for a{u) tt^tt^tt^) from [ISj. This is the major decay 
mode of the u, with B{uj — t- tt'^tt^tt^) = 0.892 ±0.007. In addition to this final state, 
we have already accounted for the 7r°7, tt+tt", and 777 final states. We thus estimate 
the contribution from oj decays not yet accounted for as being a{u — )■ 7r"'"7r~7r°) x 
S(non — Sti, H'y, r]'y)/B{u! — )■ 7r"'"7r~7r°), where ;B(non — 37r, 717, 777) = 0.0094. 

25. We have already accounted for (j) to KK, Sir, n'-f and 777. Hence there is a missing 
contribution found from i3(missing) = 1 — i3(0 — >■ KK) — B{(f) — ?■ 37r) — B{(f) — )■ 
TT'y) — B{(j) — )■ 777) = 0.0014. We then calculate a{(l)(jiorL-KK, Sn, H'y, 7/7)) = a{(p 
K+K-)B{missmg)/B{(p K+K') = 0.003 a(0 ^ K+K'). 

26. We assume that cr(7/27r^27r^) = cr(?77r^7r~27r°). 

27. One has aiKKn) = Sa{K^^K^ii^) + fr(07r°) x B{(t) KK). Here a{K^^K_^'K^) 
and cr(07r°) are obtained from BaBar (2008) In addition, B{(t) -)■ Kli) = 
0.831 ±0.003. 

28. ct{KK27i) = 9a{K+K-27i'^) + la{K+K-7i+n-). The cross sections for the i^+ir-27r 
final states are available from BaBar (2007) |19] . To estimate an uncertainty for this 
result we make use of a different procedure from [50] which uses the inclusive data 
KsX (DM1 Collaboration, [51]) as a starting point. The difference in cross sections 
between these methods is 17%, which we take to be the systematic uncertainty. 

29. We assume that a{K'^'K^Tr^n^Tr^)r^_exci = o'{K~^K~n^7v^Tr^)n-cxci- With these two 
processes as primary contributors we find a (-ft'K37r) = 2a{K~^K~7i^7r~7i'^)rj-cxd- We 
calculate a{K+K-7r+n"7r^)r,-cxci = (t{K+ K~n+n-n°) - a{(pri) x ^(0 ^ K+K-) x 
B{r] 7r+7r-7r°), where -B(0 K+K^) = 0.489 ± 0.01 and B{7] tt+tt-tt^) = 
0.2274 ± 0.0028. We obtain a{K+K^n+TT-n^) from BaBar (2007) [S3|, and (T(0r/) 
is taken from BaBar (2008) |12]. 
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When using the cross sections in the sum rules one needs the so-called undressed (or 
bare) cross-sections, ctq. The bare cross sections are obtained by removing initial-state 
radiation, but not final-state radiation, and the contribution due to vacuum polarization. 
The leptonic part of the vacuum polarization has been subtracted already from the BaBar, 
SND, and CMD2 data. However, most of the data have not been corrected for the hadronic 
vacuum polarization. We do this by using 



where n'(s) is the vacuum polarization function n'(s) = n(s) —11(0) and n(s) is obtained 
from a dispersion relation integrating over hadronic data. 
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